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ABSTRACT 

The Chandra X-ray data of the central Virgo cluster are re-examined to reveal a triple-ring structure 
around the galaxy M87, reminiscent of the spectacular triple-ring pattern of the SN1987A in the Large 
Magellanic Cloud (LMC). In the sky plane, the two apparent smaller ellipses are roughly aligned along 
the M87 jets; the larger ring centers at the M87 nucleus and is likely a circle roughly perpendicular 
to the M87 jet. Certain similarities of these two triple-ring structures might hint at similar processes 
that operate in these two systems with entirely different sizes and mass scales. We suspect that a 
major merging event of two galaxies with nuclear supermassive black holes (SMBHs) might create 
such a triple-ring structure and drove acoustic and internal gravity waves far and near. The M87 jets 
are perhaps powered by a spinning SMBH resulting from this catastrophic merging event. 
Subject headings: X-rays: galaxies: clusters — galaxies: structure — galaxies: individual (M87) — 
galaxies: clusters: individual (Virgo) 



1. INTRODUCTION 

Many similar astrophysical phenomena happen on to- 
tally different spatial and temporal scales or energy 
and mass scales. Examples include collimated jets 
from active galac tic nuclei (AGNs) containing SMBHs 
ijBegelmarJ 12003ft , fro m microquasars containing stel- 
lar mass black holes llMirab^l ^ I^odriguez? I1999D and 
from young stellar objects l|Rav et al. 1996) ; dynami- 
cal roles of waves in spiral galaxies ( Lin & Shul Il964t 
iFan fc LoiJ .199(1: Lou Fan 1998'1. in planetary rings 



lIGoldreich 



l|Fabian et al.l 



'remaind Il978() and in galaxy clusters 
2003a(l : similar high-temperature atmo- 



spher es of the Sun and of stellar mass black hole sys- 
tems (| Zhang Let al. 2000) ; similar gamma-ray flashes and 
bursts from explosio ns of perhaps m assive stars at cos- 
mological distances (jMeszarosI 120011) . from solar flares 
(^Igig£l^^^y_^9^J)^ and from the Earth's atmosphere 
l|Fishman et al.lll994 : Feng et al.lf20n2D . We report here 
the detection of an X-ray triple-ring structure in the core 
of the Virgo galaxy cluster, most likely associated with 
the AGN of the galaxy M87 (NGC 4486) and its pow- 
erful jet. This triple-ring structure is reminiscent of the 
spectacula r triple-ring pa ttern of the supernova 1987A in 
the LMC llBurrowFin 99l) . 

M87 is an active galaxy with relativistic jets in the 
central Virgo cluster. The X-ray morphology and spec- 
troscopy of the Virgo cluster have been studied pre- 



viously using data from Chandra I Young et al.l l2002() 
and XMM-Newton ijBohringer et al.ll200H iBelsole et al.l 
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20011). Comb i ning C handra. XMM and ROSAT data, 
Forman et all l)2004 hereafter F04) studied the Virgo 
cluster on various scales to identify cooling flow quench- 
ing by an AGN energy input. They noticed two prolate 
features around the M87 jet and counter jet, and referred 
to them as "cavities" caused by plasma expansions. 

In this Letter, we re-analyze the Chandra data and 
focus on the core region around M87. Two smaller ring 
structures are identified to encircle the "cavities" noticed 
by F04; a third larger ring centers at the M87 nucleus. 
We adopt the distance to M87 as 16 Mpc with 1" for 
~ 78 pc in the image (jWhitmore et al. 1995) . 

2. DATA ACQUISITION AND ANALYSIS 

We process the data of two observations (Obs-ID 2707 
and 352) pointed at M87 in the central Virgo cluster 
by the ACIS-S (Advanced CCD Imaging Spectrometer 
Spectroscopic array) instrument aboard NASA's Chan- 
dra X-ray Observatory for a total effective exposure of 
119.2 ks. The data are screened for flares where count 
rates are at least i-a away from the mean rate in the SI 
chip for 2.5-6.0 keV. The exposures before/after screen- 
ing are 98.7/89.5 ks for Obs-ID 2707 and 37.7/29.7 ks for 
Obs-ID 352. These two screened data sets are merged us- 
ing the merge-all script in CI AO 3.0.1, and further cor- 
rected with an exposure map. Another two long obser- 
vations of M87 (Obs-ID 3717 and 1808) are not included 
for analysis because of serious flare contamination. 

We avoid the over-brightness of M87 and the jet by ex- 
tracting the normalized brightness in the range [0, 8%] 
and renormalizing linearly to [0, 1], i.e., let Bi and B2 
be the brightness maps before and after the adjustment, 
we take < i3i < 0.08 from a normalized brightness 
map and renormalize it with B2 = -Bi/0.08. The 0.5- 
2.5 keV image of the central Virgo cluster is displayed in 
Figure^ by a Gaussian smoothing of 1.5" FWHM. To 
reveal structures of smaller scales in the image, we ap- 
ply an unsharp mask. A smoothed image G is obtained 
by convoluting the raw image / with a Gaussian func- 
tion. The sharpened image t/ is a weighted subtraction 
of / and G, viz. U — I -\- a{I — G) where a is a con- 
stant. By trials, we choose a Gaussian function of 15" 
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>8.0% 



TABLE 1 

Fitting parameters and orientations of the triple rings ^ 



projected ellipses 



reconstructed rings 
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ring 1 


187.7021, 12.3920 


1.9 


1.1 


76 


55 


61 


55 


ring 2 


187.7100, 12.3901 


1.5 


1.3 


70 


30 


51 


33 


ring 3 


187.7059, 12.3916 
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37 
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187.7059, 12.3911 
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a: semi major axis; b: semi minor axis; d: Angle from the major 
axis (or M87 jet direction) to the local hour circle line; a: angle 
from the ring normal (first 3 rows) or M87 jet direction (the last 
row) to the LOS; "I>: angle between the ring normal and the M87 
jet (43° to LOS); same with <I> but taken the jet orientation as 
15° to LOS. 



FWHM with a = 500%. The resulting image of Fig- 
ure is obtained by smoothing the sharpened image 
with a Gaussian function of 1.5" FWHM and by adjust- 
ing the contrast to enhance visibihty of the triple-ring 
feature; the contrast enhancement is done by displaying 
only the image within the brightness range of 15.5-24%, 
i.e., 1.24-1.92% in the original image. 

For spectra along the rings, a hardness ratio map is 
shown in Figure |21 for the central Virgo cluster where the 
hardness ratio is defined as H = (c2 — Ci)/(c2+ci) with ci 
and C2 for photon counts in the lower and higher energy 
bands, respectively (Li 2001). We select 0.5-1.0 keV and 
1.0-2.5 keV as the lower and higher bands. The images 
are smoothed with a 2" FWHM Gaussian function before 
creating a hardness ratio map in Figure [3 



3. THE TRIPLE-RING STRUCTURE 

Two nested X-ray ellipses (dashed lines marked as ring 
1 and 2 in Figure are revealed around the M87 jets. 
A third larger ellipse can be partially seen (dashed line 
marked as ring 3 in Figure ^p) . The reality of ring 1 
and 2 is evident by Figure^. Besides the northern and 
southern diffuse arcs for the initial identification of ring 
3, the ring also passes through eight bright X-ray knots 
(marked by small open arrows) that lend further support 
for its reality. The X-ray morphology of the Virgo cluster 
is fairly complex with structures on scales from ^1 — 50 
kpc (F04). We mainly focus on the triple- ring structure. 

Each presumed circular ring fits a projected ellipse us- 
ing five parameters (see TableQ, viz., the right ascension 
(RA) and declination (Dec), semi major and minor axes, 
and the inclination angle of the semi major axis to the 
local hour circle line. Because the three X-ray rings can 
not be readily recognized by our software, errors in these 
fitting parameters are estimated empirically. We note 
that position uncertainties are < 1 — 2" (~ 0.1 kpc) and 
angular uncertainties are < 5°. 

It is possible that the observed ellipses are projections 
of circular rings in the sky plane. We reconstruct 3-D ori- 
entations of the three circular rings. The jet or ientation 
may be taken as either 43° (Birctta et al.lll995l) relative 
to the line of sight (LOS) awa y from us or more likely 
15° (within 19° as inferred bv .Biretta et alJl!999) . The 
angle of the projected jet to the local hour circle line is 




1 .4% 



<1 .2% 



Fig. 1. — Two Chandra 0.5—2.5 keV X-ray images of the central 
Virgo cluster with color bars indicating the fractional flux, (a) 
— the core image plotted in linear intensity scale, smoothed by a 
Gaussian function of 1.5" FWHM and with a contrast adjustment 
to avoid the over-brightness of M87 nucleus and jets; (b) — the 
triple-ring feature (dashes) enhanced image with an unsharp mask. 
In the sharpened image, two rings (northwest ring 1 and southeast 
ring 2) are roughly along the M87 jet. The third ring 3 can be seen 
in broken segments, viz. the northern and southern portions fit in 
an ellipse with the M87 nucleus as its center. Ring 3 passes across 
eight bright X-ray knots identified by small open arrows. 



~ 70° (FigureP). Thus, the angle $ (see TableUJ of the 
reconstructed circular ring normal to the jet direction (if 
43° to LOS) is 61°, 51° and 9° for ring 1, 2 and 3, re- 
spectively. This means that ring 3 is likely a circular ring 
nearly perpendicular to the M87 jet. For a jet orientation 
of 15° to the LOS, the corresponding angles of ring 
normal to the jet (see Table are 55°, 33° and 22°, re- 
spectively. The derived values of $ and for presumed 
'circular' ring 1 and 2 are fairly large; these two rings 
might be elliptical in the 3-D space or their orientations 
are not perpendicular to the jet. Our estimate excludes 
the possibility that the elliptical shape of rings 1 and 2 
is caused by relativistic projection effects if the two cir- 
cular rings expand radially away from the jet axis while 
moving apart from each other relativistically along the 
line connecting the two ring centers, because by special 
relativity the ring would be elongated transverse to the 
jet direction. We note that ring 3 passes through eight 
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X-ray bright knots with its center located exactly at the 
M87 nucleus. The chance of such a coincidence is small. 

As projected ring 1 and 2 roughly align along the di- 
rection of the M87 jet and counter jet, their intrinsic 
elliptical (instead of circular) shapes in 3-D space are 
plausible. Given relatively small angle of the ring 3 nor- 
mal to the jet (9° or 22°), one may qualitatively regard 
ring 3 as circular and perpendicular to the jet. 

4. SPECTROSCOPY OF THE THREE RINGS 

From Figure 13 the softest regions locate at (i) around 
the jet within a ~20" radius sector from west to north 
of the nucleus, (ii) an extended arc coincident with the 
east portion of ring 1 from north to south of the nu- 
cleus, (iii) a short arc coincident with ring 3 about ~30" 
away from the nucleus in the northwest and (iv) some 
extended filaments southeast of the nucleus outside the 
three rings. By comparing Figure [5] wit h the Ha-|-[N 
h] em issions of the central Virgo cluster ijSoarks et al.l 
I1993D . we find that the soft region (i) and (iv) are corre- 
lated with Ha-|-[N h] filaments. Si milar correlations wer e 
also found in the Perseus cluster ijFabian et alJ 12003131) . 
indicating that intracluster gases become cooler when ad- 
jacent to the filaments due to thermal conduction. 

The association of regions (ii) and (iii) with ring 1 and 
3 respectively is interesting with region (ii) being more 
striking. A harder spectrum may correspond to a higher 
temperature. Different temperatures along ring 1 indi- 
cate that X-ray emissions may not be simply caused by 
a high density of electrons alone, but may involve super- 
sonic shock flows from the AGN. Other parts of rings 
without temperature difference relative to surroundings 
might indicate subsonic flows. 

Energy spectra are derived from two segments of ring 
1, viz., the east part identified with region (ii) and the 
west part. To compare the spectra in these two parts 
with their environments, the source region is partitioned 
as an elliptical annulus along ring 1 with a width of 
and a background region of a bigger elliptical annulus 
outside 5" of ring 1 with a width of ~2". The source 
and background regions are separated into two parts by 
a straight solid line in Figure |21 

In Figure |21 the east part carries features, while the 
west part blends with the environment. We adopt a hot 
diffuse gas emission model (VMEKAL) to fit the spectra 
of both parts, with free parameters of equivalent hydro- 
gen column, temperature and elemental abundances of 
O, Ne, Mg, Si, S and Fe. The fitting is carried out using 
XSPEC 11.3. The best-fit temperatures are 0.74t°:EJ| 
keV (xVd.o.f. 71/54) for the east part and 1.321°:°^ 
keV (x^/d.o.f. = 53/55) for the west part. A remarkable 
difference in plasma temperature can be seen from energy 
spectra using the hardness ratio map in Figure|21 In con- 
trast, only the central ring of SN1987A c an be resolved 
by Chandra using a sub-pixel technique (Bu rrows et al.l 
iOOO; Park ct al. 2002; Michael ct al. 2002). As its X-ray 
spectrum may be fit by a plane-parallel shock model, the 
central ring of SN1987A is thought to involve an expand- 
ing blast wave. By this clue, we attempt to fit the spectra 
of both east and west parts of ring 1 (FigureOJ)) using the 
plane-parallel shock model (VPSHOCK) with the same 
free parameters of elemental abundance indicated above. 
The west part spectrum does not fit the plane-parallel 
model well with a large x^/d.o.f. ~ 5, while the east part 



3 




Fig. 2. — Hardness ratio map in the central Virgo cluster. The 
two energy bands are 0.5-1.0 keV and 1.0—2.5 keV. Each pixel 
stands for H = (c2— ci)/ (c2-|-ci) where ci and C2 are photon counts 
at this pixel in the lower and higher energy bands respectively. Four 
softest regions are labeled. The solid line indicates a separation of 
the two parts along ring 1 for spectral comparisons (see §|4j. 



spectrum fits well with the shock model. The best-fit 
temperature is 0.7lto;o4 keV (x^/d.o.f. = 77/53), consis- 
tent with the hot plasma model. This suggests that ring 1 
may involve a blast wave where the west part might have 
been thermalized with the surroundings and the east part 
might be in the process of being thermalized. For the east 
part, the best-flt abundances of elements O, Si, S, and Fe 
in the thermal model are systematically lower than that 
in the west part. The enhancement of abundances in 
the warmer region may indicate more active nuclear pro- 
cesses in the past, perhaps caused by the thermalization 
in the west part. The similar temperature-abundance 
correlation is al so found with XMM in the inner region 
for radius <1' (IBohri nger et al J 1200 ID and for the two 
large-scale arms ( Bclso lemLLll200 1(). Because of uncer- 
tainties in estimating abundances with Chandra data for 
the fine structures, we refrain from more detailed com- 
parisons with previous results of XMM for the larger 
structures. 

5. DISCUSSIONS 

The triple-ring pattern revealed in the central Virgo 
cluster is the only one obs erved since the discovery of 
triple rings in SN1987A8 Ib urrow ^ 11991 . The ring 
sizes in the Virgo cluster are several thousand times 
larger than those of SN1987A. The triple-ring structure 
of SN1987A has been simulated in terms of interact- 
in g winds near the end of st ellar evolution. According 
to lTanaka fc Washimil l)2002D . during the star's red su- 
pergiant phase, slower winds persist. Prior to the su- 
pernova explosion, the red supergiant somehow evolved 
into a blue supergiant that drove faster winds. Mag- 
netohydrodynamic (MHD) interactions of fast and slow 
winds from a rotating star produce the triple-ring struc- 
ture in SN1987A. We suspect that a triple-ring structure 
might be generic in a rotating system involving magne- 
tized faster winds catching up slower winds (Lou l|)9j, 
Il996j) . Meanwhile a central source of radiation is required 

* http:/ /opositc.stsci.cdu/pubinfo/jpcg/SN1987A_Rings.jpg 
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for the rings to shine in the X-ray band. For SN1987A, 
the strong UV radiation from the supernova provides the 
necessary illumination. 

We suspect that a catastrophic merging event around 
the M87 nucleus might be responsible for the triple-ring 
structure revealed here. Merging of two SMBH systems 
has been proposed to drive the o bserved X-shaped ra- 
dio lobes l)Merritt fc Ekersl |2002|) . In this scenario, a 
'slower wind' was present before the merging begins, e.g., 
a mixture of galactic winds from two merging galaxies. A 
'faster wind' was then driven during the merging process. 
It has been suggested that an active AGN phase o ccurred 
in th e core region of M87 at about 10^ years ago ijKaiseil 
|2003|) : perhaps, this AGN was triggered by the merging 
event discussed above, when significantly higher accre- 
tion rate was provided by the merger. The merger may 
fuel both the radiation energy output of the AGN and 
the growth of the supermassive black hole in the center; 
huge amount of high energy radiation shines the triple 
rings in X-ray bands. It is also plausible that the result- 
ing SMBH s pins rapidly to pow er the highly coUimated 
M87jets ( Blandford fc Znaiekl [r9771. 

F04 revealed two expanding cavities around the M87 
jet and counter jet, enclosed within our ring 1 and 2 re- 
spectively. By different external pressures, they inferred 
the far front of the cavities is mildly supersonic while 
the inner part, close to the nucleus, is collapsing. We 
suggest that these two rings are actual rings instead of 
projections of prolate spheres, because when two cavities 
collide into each other it is hard to maintain their spheri- 
cal shapes in the spatially overlapping portions of the two 
spheroids, and consequently their projections cannot be 
viewed as two complete and nested rings. We further sug- 
gest that when the jets form, they push materials around 
and drive powerful blast waves of ring shapes. 

Qualitatively similar central str uctures are present 
in the clusters Perseus llPabian et al.l l2008aD. Hydra 
A llMcNamara et al.ll200lD Centaurus A (jKraft et all 
120021 iKarovska et all 12002) and Virgo (F04). Physi- 
cally, core activities can excite large-scale acoustic waves 
(p— modes) and internal gravity waves (5— modes) in a 
deep er avitational potential well of central galaxy cluster 
ifiTiTiiri 99,^^1 . Cluster core activities can be either violent 
and sporadic or sustained at a level or a combination of 
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both. Depending on the phase of acoustic wave propa- 
gation and on physical properties of core environment, 
we may observe various X-ray morphologies superposed 
onto a smooth luminous core. Beautiful a coustic rings or 
arcs have been revealed in Perseus (Jabian..et al. 2003^ 
and in Virgo on larger scales (F04). Intense core X-ray 
emissions tend to induce inflows of gas that perturb the 
dark matter via gravity. The violent relaxation jaroceeds 
in a sound-crossing timescale |Lvndcn-Bcll 1967). This 
could be a persistent source of modes and g— modes in 
a cluster core. By such wave and Laudau-damping pro- 
cesses, slowly inward drifting gas and dark matter are 
virialized to sustain X-ray losses from hot electron gas. 
As a cluster core may involve magnetic field of strengths 
up to ^ 30 — 40/xG, MHD waves should also play sig- 
nificant roles. Based on the hardness ratio map and 
the energy spectra fitting, we infer only parts of ring 
1 and 3 involve supersonic shock flows and are in pro- 
cesses of being thermalized with the surroundings. By 
this reasoning, other parts without apparent difference 
with surroundings might have been thermalized already. 
In summary, we hope that the triple-ring structure of 
M87 in the central Virgo cluster would stimulate more 
theoretical, simulation and observational studies of MHD 
flow, wave and shock interactions. 
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